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Anderson, et al. find the measured trajectories of Pioneer 10 and 11 spacecraft deviate from the 
trajectories computed from known forces acting on them. This unmodelled acceleration can be 
accounted for by non-isotropic radiation of spacecraft heat. Various forms of non-isotropic radiation 
were proposed by Katz, Murphy, and Scheffer, but Anderson, et al. felt that none of these could 
explain the observed effect. This paper calculates the known effects in more detail and considers 
new sources of radiation, all based on spacecraft construction. These effects are then modelled over 
the duration of the experiment. The model provides a reasonable fit to the acceleration from its 
appearance at a heliocentric distance of 5 AU to the last measurement at 71 AU, but overpredicts 
by 9% the decrease in acceleration between intervals I and III of the Pioneer 10 observations. (For 
comparison, the two different measurements of the effect (SIGMA and CHASMP) themselves differ 
by 4% in interval III.) In any case, by accounting for the bulk of the acceleration, the proposed 
mechanism makes it much more likely that the entire effect can be explained without the need for 
new physics. 
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I. INTRODUCTION 



II. THE ANOMALOUS ACCELERATION 



In Anderson et al. compare the measured trajec- 
tory of several spacecraft against the theoretical trajec- 
tory computed from known forces. The find a small but 
significant discrepancy, referred to as the unmodelled or 
anomalous acceleration. It has an approximate magni- 
tude of 8 X 10~^ cm s~^ directed approximately towards 
the Sun. Needless to say, any acceleration of any object 
that cannot be explained by conventional physics is of 
considerable interest. Explanations for this acceleration 
fall into two general categories - either new physics is 
needed or some conventional force has been overlooked. 

One of the most likely candidates for the anomalous 
acceleration is non-isotropic radiation of spacecraft heat. 
This is an appealing explanation since the spacecraft dis- 
sipates about 2000 watts total; if only 58 watts of this 
total power was directed away from the sun it could ac- 
count for the acceleration. Several possible mechanisms 
have been debated in the literature, but none are totally 
satisfactory. 

In this paper we re-examine each proposed mechanism, 
explicitly including their time dependence. We propose 
several additional mechanisms - asymmetric RHU heat, 
misdirected feed radiation, and mis-modelled solar reflec- 
tivity. Finally, we compare the acceleration induced by 
the proposed mechanisms with the measured data, and 
get reasonable agreement over the whole data span. 



'Electronic address: lou@cadence.com 



As the Pioneer spacecraft receded from the sun, solar 
forces decreased and only gravitational forces, and an oc- 
casional maneuver, affected the trajectory of the space- 
craft. Anderson, et al. noticed that a small additional 
acceleration needed to be added to the known forces to 
make the measured data and computations match. This 
is the anomalous acceleration, which started to become 
noticeable about 5 AU from the sun, and was roughly 
the same for Pioneer 10 and 11. The onset is shown in 
Figure [l[ 

Further constraints come from the ongoing study of Pi- 
oneer 10, where there are fewer confounding effects and 
the data span is long enough to provide significant con- 
straints due to the radioactive decay of the heat sources. 
Figure |^, reproduced from j2| , shows the measured ac- 
celeration 1987 to 1998. (Although they have different 
horizontal axes. Figure ^ largely follows Figure |l| chrono- 
logically. Pioneer 10 was at 40 AU in 1987.) The authors 
divide the 1987-1998 Pioneer 10 history into three inter- 
vals. Interval I is January 1987 to July 1990, interval II 
from July 1990 to July 1992, and interval III is from July 
1992 to June 1998. The authors make this distinction by 
looking at the spin rate of the craft - in intervals 1 and 111 
it was decreasing smoothly, but in interval II it decreased 
quickly and irregularly. They therefore consider the data 
from interval II to be less reliable than intervals I and 
III, since whatever affected the spin (probably gas leaks) 
may also have affected the acceleration. 

More recent analyses have refined these results some- 
what, though the main conclusions remain unchanged. 
Table | shows the most recent results from which 
fits a constant, independent acceleration in each interval. 
Accelerations are in units of 10~* cm s~^. SIGMA and 
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CHASMP are two different and largely independent tra- 
jectory modelling programs; the difference between the 
programs is our best estimate of the real uncertainties 
since it is far greater than the formal errors. This data, 
taken at face value, shows that 57 directed watts can ac- 
count for the acceleration in 1998, and that a 3% decrease 
was observed between interval I and interval III. 



TABLE I: WeiEhted Least Squares (WLS) results from An- 
derson, et al.M, and equivalent directed power for a 241 kg 
spacecraft mass 



Interval 


SIGMA 
accel. 


watts. 


CHASMP 
accel. 


watts 


Jan 87- Jul 90 


8.00 ± 0.01 


57.8 


7.84 ± 0.01 


56.7 


Jul 92 - Jul 98 


8.25 ± 0.03 


59.6 


7.91 ± 0.01 


57.2 



III. PREVIOUS WORK 

Many paper |Q and web||] descriptions of the Pioneer 
spacecraft are available. In this section we summarize the 
existing literature on the hypothesis that non-isotropic 
radiation is responsible for the unmodelled acceleration. 

Murphy 1^ (and a related proposal by Scheffer[0) sug- 
gests that the anomalous acceleration seen in the Pi- 
oneer 10/11 spacecraft can be, "explained, at least in 
part, by non-isotropic radiative cooling of the space- 
craft." KatzQ proposes that at least part of the acceler- 
ation is generated by radiation from the RTGs reflecting 
off the back of the antenna. Slusher (as credited by An- 
derson) proposed that the forward and backward surfaces 
of the RTGs may emit non-equally. Anderson, et al. ar- 
gue in reply |, |, that none of these proposed sources 
adequately account for the acceleration. 



A. Radiation of spacecraft power 

First, consider thermal radiation from the body of the 
spacecraft. A thought experiment shows that the elec- 
trical power dissipated in the spacecraft must result in 
thrust. The simplest model consists of the main com- 
partment as a 60 watt isotropic radiator, and the back 
of the antenna a mirror. The antenna subtends 120 de- 
grees as seen from the instrument compartment, so if 
the emitted radiation is isotropic, the antenna intercepts 
1/4 of the total radiation, and reflects it away from the 
sun. Since the main compartment is centered behind the 
antenna, and since the sides, if anything, are worse ra- 
diators than the front, we conclude that at least 25% of 
spacecraft electrical power must be converted to thrust. 

A more detailed analysis shows the radiation is even 
more anisotropic than these arguments would suggest. 
Assuming a uniform internal temperature and closed lou- 
vers, the power emitted from each surface is propor- 
tional to the area times the "effective" emissivity of the 
surface |l^. The sides and the rear of the compartment 
are covered with multi-layer insulation(MLI)|Q, with an 
effective emissivity of 0.007 to 0.01 |Q. The lowest emis- 
sivity material on the front of the spacecraft is the sur- 
face of the louvers, with a emissivity of 0.04Q. Since the 
sides and the front have comparable surface areas, then 
about 80% of the total power will be radiated though 
front (though hard to characterize heat leaks could re- 
duce this value). Frontal radiation would be expected to 
be about 66% efficient, assuming Lambertian emission. 
Defining eBus as the fraction of main compartment heat 
that is converted to thrust, we then expect esus to range 
between 0.25 (blockage arguments) to 0.52 (differential 
emissivity) . 

From , the total electrical power is modelled 
E{d) ^ (68 + 2.6 (1998.5 - d)) watts 
and the thrust (assuming an 8 watt radio beam) is 
BUS{d) = EBUS {E{d) - 8.0 watts) 



IV. DISCUSSION 

We consider asymmetrical radiation from 4 sources - 
the RTG heat (direct radiation and reflection off the an- 
tenna), the electrical power dissipated by the spacecraft, 
the radioisotope heater units (RHUs) on the spacecraft, 
and radiation from the feed that misses the antenna. We 
also consider one modelling error, a mis-estimation of the 
reflectivity of the antenna to solar radiation. The avail- 
able power from all these sources changes in time. In 
the following discussion, let d be the date, in years. The 
sunward side of the spacecraft is the back, and the anti- 
sunward side, in the direction of motion, is the front [pl]| . 
We calculate thrust in units of watts of directed (anti- 
sunward) radiation. 



B. Feed pattern of the radio beam 

An ideal radio feed antenna would illuminate its dish 
uniformly, with no wasted energy missing the dish. How- 
ever, the feed is physically small and cannot create such a 
sharp edged distribution, so some radiation always spills 
over the edge. This radiation is converted to thrust with 
an efficiency of 1.7 since it directly subtracts from the 
sun directed power and adds anti-sun power at a roughly 
45 degree angle to the spin axis. This produces thrust 

RADIO{d) = epEED (8 watts) 1.7 

where cfeed is the fraction of RF power that misses the 
antenna. Since dish area is wasted if not fully illumi- 
nated, an optimum feed (for transmission) will result in 
^FEED ~ 0.1. 
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C. Radiation from the RHUs 

From diagram 3.8-1 in |Q, 10 1-watt (in 1972) radioiso- 
tope heater units are mounted to to external components 
(thrusters and the sun sensor) to keep them sufficiently 
warm. The diagram is not very specific, but the units to 
which they are mounted are primarily behind the main 
dish. Radiation from these components will contribute 
thrust, which we model as 

RHU{d) = €RHu (10.0 watts)2-('^-i^^^)/88 

where €rhu is the proportion of RHU heat converted to 
thrust. Reasonable values for €rhu might range from 
0.0 to 0.5, with the latter corresponding to components 
behind the dish radiating uniformly. 

D. Radiation from the RTGs 

The RTGs might contribute to the acceleration by 
radiating more to the front of the spacecraft than the 
rear, and/or by having their heat refiected asymmetri- 
cally from the spacecraft. The RTGs radiate all the ther- 
mal power that is not turned into electricity, so 

RTGMEAT{d) = (2580 watts)2"'''~i^^2)/^^ - E{d) 

In direct radiation asymmetry is estimated to con- 
tribute to thrust with an efficiency of at most ±0.003. 
RTG reflection by the antenna was proposed by Katz, but 
argued against by Anderson, primarily on the grounds 
that the RTGs are on-axis as seen by the antenna. We 
re-examine this argument here. From figure 3.1-2 of 
we see that the centerline of the RTGs is behind the 
center of the antenna. Measurements from this diagram 
indicate this distance is about 23.8 cm. Figure 3.1-3 of 
g shows the far end of the RTGs is 120.5 inches (or 3.06 
meters) from the centerline. From this geometrical data 
we can estimate the area blocked by the antenna from 
each RTG[|lJ]. Numerical integration of these areas, as- 
suming Lambertian emission by the RTGs, shows about 
0.6% of the near RTG radiation and 0.4% of the far RTG 
radiation fall upon the dish. This energy is turned into 
thrust by two effects. First, the antenna shadows radia- 
tion which would otherwise go forward. An angle in the 
middle of the antenna is about 17 degrees forward; this 
corresponds to an efficiency of 0.3 (the true efficiency is 
probably higher since the edge is both at a greater angle 
and more brightly illuminated.) Next, the energy that 
hits the antenna must go somewhere. Some will be ab- 
sorbed and re-radiated; some will bounce into space, and 
some will bounce and hit the instrument compartment, 
and be reflected or re-radiated from there. A detailed ac- 
counting seems difficult, but an overall efficiency of 0.7- 
0.9 seems reasonable (0.3 for shadowing and 0.4-0.6 for 
refiection and re-emission). 

We model the total thrust from RTG heat as 

RTG{d) = ERTG RTGMEAT{d) 



where tRTC is the proportion of RTG heat converted to 
thrust. Gombining the effects of this section, we expect 
(-RTG to range from 0.004 to 0.012. 



E. Antenna solar reflectivity 

The trajectory analysis programs fit the reflectivity of 
the spacecraft to solar radiation, /C, as a force that falls 
off as l/r^, where r is the heliocentric distance. This flt 
can hide an otherwise unmodelled acceleration. Over a 
short time period, during which r varies little, any con- 
stant radial acceleration can be absorbed into /C. Over 
a longer period of time the fitting procedure will mask 
any component of anomalous acceleration that varies as 
1/r^ and is less than the acceleration corresponding to 
the allowed variation in /C. In particular the acceleration 
proposed in this paper will be partially masked since it 
decreases with time and hence has a 1/r^ component. 

The fitted solar reflectivity constant also provides a 
natural explanation for the onset of the anomalous accel- 
eration. Consider the case where the acceleration (from 
any cause) exists for all r. When r is small, the fitting 
programs absorb the extra acceleration by adjusting the 
value of /C. As r increases, the power available from this 
source decreases, and eventually IC runs into the limits 
allowed in the flt. (Physically reasonable values perhaps 
range from 1.5 to about 1.9; they are certainly greater 
than 1.0 and less than 2.0). Once the limit of adjust- 
ment for K, is reached, it becomes constant and can no 
longer mask the acceleration, which appears as shown in 
flgure ||. It might be possible to see additional signs of 
this process in archival data - it would show up as a de- 
crease in the fitted value of IC as the spacecraft receded 
from the sun. 

In this paper, we model the effect of any error in K. 
by introducing a flctitious force, whose value is simply 
the solar force on the spacecraft times the error in IC. 
We assume the distance from the sun, measured in AU, 
increases linearly from 20 AU in 1980 to 78.5 AU in 2001: 

r{d) = 20 -f (d - 1980)/21 • (78.5 - 20) 

The thrust, in watts, is 

SOLAR{d) = KsoLAR ^(1.37 Ynffe/r^^d) 

where = 1367 W/m'^(AU)^ is the "solar radiation 
constant" at 1 AU, and Ksolar, the amount by which 
the solar reflection constant is underestimated. 



V. COMPARISON WITH EXPERIMENT 

To compare the hypothesis with experiment, we sum 
the individual sources, then convert to acceleration by 
dividing by c, the speed of light, and m, the spacecraft 
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mass (here 241 kg): 



50-Day Acceleration Averages for Pioneer 1 

Dashed Cur\/e is Best— Fit Anomalous Acoeleration 



acc{d) 



1 



■[RHU{d) + RTG{d) 



c • m 

RADIO{d) + BUS{d) - SOLAR{d)] 



We then compare with the plots from |j . The pro- 
posed explanation has 5 adjustable parameters. In the- 
ory all are separable since they decay at different rates; 
in practice the data are not good enough to separate 
them and many fits are plausible. One reasonable fit 
over the entire data span has the following coefficients: 
<^RHU = 0.5, eRTG = 0.0108, epEED = 0.1, esus = 0.35, 
and KsoLAR = 0.3. 

This fit to the data is shown in Figures |l|and|. The 
agreement seems reasonable in both regimes, and the pro- 
posed model provides a better fit to the early data than 
the constant acceleration of ^ , even assuming reflectivity 
mismodelling to account for the onset of the acceleration. 
The fit from 1987 to 1998 also looks acceptable, as shown 
in Figure |^. 




FIG. 1: Anomalous acceleration vs heliocentric distance in 
AU, from (error bars), model prediction from this paper 
(solid line), and constant acceleration plus reflectivity mod- 
elling error(dashed line) 



Finally, we compare with the most recent results ||] 
that fit a constant acceleration in each interval of the 
later Pioneer 10 data. The proposed model gives an av- 
erage thrust of 57.8 watts in interval I, and 51.0 watts in 
interval IIL We can normalize the result to get the cor- 
rect overall average, or the right acceleration in interval 
I, but in either case we would expect to see an 11.8% 
decrease from interval I to III, where only a 3% decrease 
is observed. The two different measurements of the effect 
(SIGMA and CHASMP) themselves differ by 4% in in- 
terval III. If we treat this difference as a statistical result 
(a procedure of dubious merit, but the best we can do) 
then the 9% discrepancy is 2.25 standard deviations out. 
This makes it unlikely at about the 2% level that this 
hypothesis alone accounts for all the measured result. 

We can get a better fit (1.75 sigma) to the Pioneer 
10 data by assigning different efhciencies to instrument 
heat and main compartment heat, at the cost of an extra 




1987 1938 1983 1990 1991 1992 19S3 1994 1995 1996 1937 1993 1999 

Date 

FIG. 2: Figure from with fitted data added. The dotted 
line is Turyshev's empirical fit; the solid line is the model 
hypothesized in this paper. 



parameter and the need to consider instrument power 
dissipation in detail [H. 



VI. CONCLUSIONS AND FUTURE WORK 

There is surely an unmodelled effect on the Pio- 
neer spacecraft, based upon its thermal characteristics. 
Rough estimates show it can account for the magnitude 
of the unmodelled acceleration to within the errors, but 
overpredicts the rate of change. In any case, the pro- 
posed explanation, by accounting for the bulk of the ef- 
fect, makes it more likely that conventional physics can 
account for the entire unmodelled acceleration. Conven- 
tional explanations for the remaining discrepancy include 
other unmodelled effects such as gas leaks, inaccuracies 
in the simple thermal model, or the effects of a complex 
fitting procedure applied to noisy data. 

This explanation also explains some other puzzles: the 
values of acceleration of Pioneer 10 and 11 would be ex- 
pected to be similar, but not identical, as observed. The 
acceleration would not have a strong effect on the spin; 
most of the radiation will generate little torque. Other 
spacecraft, built along the same general principles, would 
be expected to show a similar effect, but planets and 
other large bodies would not, as is observed. 

More detailed modeling, using the Pioneer materials, 
construction details, and history, might confirm or refute 
the proposed hypothesis, and additional tracking could 
be useful as well. However, such improvements are lim- 
ited since accurate thermal modelling is difficultly and 
the spacecraft was not designed for this purpose. Longer 
term, other proposed experiments such as LISA||l^ are 
designed specifically to reduce non-gravitational system- 
atics (by a factor of about 10^) and allow frequent and 
accurate tracking (a differential distance measurement. 
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each second, accurate to 10~^ cm) Assuming the anoma- 
lous acceleration exists at all heliocentric distances, (as 
argued in section IV E), then it should be detectable in 
just a few seconds of LISA data. On the other hand, if no 
unmodelled acceleration is detected in these more precise 
experiments, then almost surely the anomalous acceler- 
ation of Pioneer 10/11 is caused by overlooked prosaic 
sources such as those proposed here. 
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